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Abstract 

We estimate the effects of viscosity on the phase space distribution appearing in the Cooper-Frye 
formula within the framework of the Grad's fourteen moment method and find that there are non- 
trivialities in the discussion of a multi-component system. We calculate the viscous corrections 
of particle spectra and elliptic flow coefficients from the distortion of the distribution using the 
flow and the hypersurface taken from a (3+l)-dimensional ideal hydrodynamic simulation. We 
see that the bulk viscosity have visible effects on particle spectra. The results suggest we should 
treat both shear and bulk viscosity carefully when constraining the transport coefficients from 
experimental data. 



1. Introduction 

Ideal hydrodynamic models have been successful U 0] for the quark-gluon plasma (QGP) 
created at the Relativistic Heavy Ion Collider (RHIC) in Brookhaven National Laboratory (BNL). 
Our next step is to develop viscous hydrodynamic models because the success suggests the devi- 
ation from the equilibrium is not so large. We mainly consider the effects of bulk viscosity since 
the importance of bulk viscosity in terms of the violation of scale invariance near the crossover 
temperature T c has been re-recognized Q3D - Particle spectra are affected by viscosity in two ways: 
one is the variation of the flow Sif — - and the hypersurface 6dcr M - dcr^ - dcr^, the other 
is the distortion of the distribution 6f — f - fo because one needs to employ the Cooper-Frye 
formula 

M=(kL'^" if ° +sfi) m 

at freezeout to convert the hydrodynamic results into the particle picture. Here the index i denotes 
particle species and gj the degeneracy. This enables us to follow further dynamical evolution of 
the hot QCD matter in the cascade models. In this study we concentrate on corrections to particle 
spectra due to the modification of the distribution at freezeout in the framework of Grad's method. 



2. Relativistic Kinetic Theory 

We express the distortion of the distribution 6f in terms of macroscopic variables for a 
relativistic multi-component gas within the framework of the Grad's 14-moment method. For 
this we generalize the method developed in Ref. [4] for a single component system by specifying 
non-trivialities proper to a system of multi-species particles. The details are shown in Ref. ||5|]. 
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We have 14 equations that relate the macroscopic variables with the microscopic distri- 
bution 6ff. the kinetic definitions of the bulk pressure II = —^A^dT^, the energy current 
W" = A^ST^Up, the charge current V = A^SNg and the shear stress tensor = 5T^ v) = 
[jiA^A^ + A^A v a ) - ±A" V A a/j ]6T afi and the two Landau matching conditions u^ST^Uy = and 
UfjdNfi = 0. Here 6T MV and 6Ng are the dissipative parts of the energy-momentum tensor and the 
baryon number flux, is the flow and A /<v = g>™ - u^u v is the projection operator where the 
metric g^ v = diag(+, -,-,-). We need the Landau matching conditions to ensure the thermo- 
dynamic stability to the first order. This can be shown by explicitly calculating dis^u^) /dH\n=o, 
where s>* is the entropy current. These equations can be written in terms of 6f by using the 
expressions of 6T MV and SNg in the relativistic kinetic theory for a multi-component system, 



(2) 



where is the baryon number. 

The Grad's 14-moment method allows us to write the distortion of the distribution for the z'-th 
particle species with 14 unknowns as 



(3) 



where the positive sign is for bosons and the negative one for fermions. In the conventional 
formalism one has the scalar term e, but here we do not have it because the traceless condition 
is not imposed on the tensor term s^. These two formalisms are not equivalent for a multi- 
component system because the scalar term in the conventional method has no particle species 
dependence whereas the trace part in our method does in the sense that it is dependent on the 
mass of the ;'-th particle species. 

We obtain three independent sets of equations by inserting 8f into the conditions: 
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where J mn is a coefficient of the expansion of the moment 
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J mn is defined in the same way but with the baryon number bi as an additional weight factor in the 
summation over the particle species index i. The unique form of 6f can be obtained by solving 
them. The correction tensors are 



£jU = Don^+DjW^ + DiV^, 



V = (BoAfn, + BoUf,Uy)n + WiUfrWy) + 2BiU(f,V Y ) + B 2 jr fn „ 
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(9) 
(10) 



where the prefactors D, and B,- are functions of J„ m and J mn . One cannot take the scalar term in 
(0 because that induces a change of the sign in the bulk pressure-related prefactors Do, Bo and 
Bo which causes divergence in the case of the 16-component hadron resonance gas discussed 
later. Note here that (i) information of all the components in the gas comes into the viscous 
modification of a particle species through the transport coefficients, J mn 's and 7,,,,,'s, and (ii) if 
one takes the quadratic ansatz in Ref. [6], i.e., = Cin MV + C2A /JV II, the system becomes 
thermodynamically unstable because the matching conditions are not satisfied. 

3. Results 




Figure 1 : Effects of bulk viscosity on (a) p^-spectra and (b) elliptic flow coefficients V2(pt) for the negative pion. Solid, 
dotted, and dashed lines are, respectively, the results without any viscous corrections, with the effect of bulk viscosity 
with a = 15/2 and with a = 15. The thick(thin) lines show that the absolute value of the ratio of the correction to the 
ideal spectrum is smaller(larger) than 0.5. 

One needs models for the equation of state (EoS), the macroscopic variables and profiles of 
the flow and the hypersurface to estimate the effects of 6f on particle spectra. We consider 16- 
component hadron resonance gas for the EoS, which has mesons and baryons with the mass up 
to A(1232). The bulk pressure is estimated with the Navier-Stokes limit, i.e., II = -^d^ where 
the models of bulk and shear viscosity are ( — a(c 2 s - ^) 2 t] J3,§,@] and 7 = r- ifioll . Here s and c s 
are the entropy density and the sound velocity, respectively, a is a free parameter. The freezeout 
temperature T f is set to 0.160 GeV where 77 = 1.31 x 10 3 GeV 3 and £ = 4.37 x 10~ 4 GeV 3 when 
a = 15. The flow and the hypersurface are taken from a (3+l)-dimensional ideal hydrodynamic 
simulation [ 1 , 2] in Au+Au collisions at y/s^N = 200 GeV where the impact parameter is 7.2 fm. 
It should be noted again that what we estimate here is not the distortion of the flow du^ but the 
viscous corrections through the modification of the distribution Sf at freezeout. For numerical 
estimations we employ the zero net baryon density limit, which is well justified for the hot QCD 
matter at RHIC. We do not have to worry about the loss of the Landau matching condition of 
the baryon number flux because the condition still yields a finite relation even in the limit of 
vanishing chemical potential 

We estimate the particle spectra of negative pion with mass m = 0.139 GeV using the Cooper- 
Frye formula ([T]) without changes of flow velocity nor freezeout hypersurface from the ideal 
hydrodynamic results. The bulk viscosity lowers the mean pj of the particle spectra, (pr), 
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when estimated with the ideal hydrodynamic flow and hypersurface. This can be attributed to 
the reduction of the effective pressure due to the bulk viscosity. The isotropic nature of the 
bulk pressure reduces the pressure in transverse directions when the system is expanding in 
longitudinal direction. On the other hand, the bulk viscous correction enhances the elliptic flow 
coefficients V2(pr)- This can be interpreted in the relation dv2(pr)/dpT — V2/{pr) 111], because 
the slope of differential V2(pr) is expected to become steep when the (pr) is decreased while the 
average V2 is not affected by viscosity with in our calculations. 

We have to be aware that the viscous corrections might have been overestimated because (i) 
we take the Navier-Stokes limit for fl and thus do not consider relaxation effects, and (ii) the 
derivatives of ideal hydrodynamic flow is generally larger than those of real viscous hydrody- 
namic flow as viscosity tends to eliminate differences in the flow. 

4. Conclusions 

We obtained the unique and consistent expression of the distortion of the distribution for a 
multi-component system in the framework of the Grad's 14-moment method. We find several 
non-trivial issues along the course of discussion. First, one needs to introduce the non-zero 
trace tensor term in the correction instead of the conventional scalar term in the case of the 
16-component hadron resonance gas to get a finite result. Second, the distortion is different 
by particle species, but is dependent on all the components in the gas because the transport 
coefficients and macroscopic quantities J m „ and J,„„ are. Third, the matching for the baryon 
number flux remains meaningful even in the limit of vanishing chemical potential. 

Effects of bulk viscosity on particle spectra through the distortion of the distribution is es- 
timated, and we find it suppresses /^-spectra and enhances V2(pr) when estimated with (3+1)- 
dimensional ideal hydrodynamic flow and hypersurface. We should be careful of the possible 
overestimation in the usage of the Navier-Stokes limit and the ideal hydrodynamic flow. A 
(3+l)-dimensional viscous hydrodynamic model is necessary to accurately estimate the viscous 
effects to compare with experimental data and to constrain the transport coefficients from them. 
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